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We construct an exactly soluble spin-| model on a honeycomb lattice, which is a generalization of 
Kitaev model. The topological phases of the system are analyzed by study of the ground state sector 
of this model, the vortex-free states. Basically, there are two phases, A phase and B phase. The 
behaviors of both A and B phases may be studied by mapping the ground state sector into a general 
p-wave paired states of spinless fermions with tunable pairing parameters on a square lattice. In this 
p-wave paired state theory, the A phase is shown to be the strong paired phase, an insulating phase. 
The B phase may be either gapped or gapless determined by the generalized inversion symmetry is 
broken or not. The gapped B is the weak pairing phase described by either the Moore- Read Pfaffian 
state of the spinless fermions or anti-Pfaffian state of holes depending on the sign of the next nearest 
neighbor hopping amplitude. A phase transition between Pfaffian and anti-Pfaffian states are found 
in the gapped B phase. Furthermore, we show that there is a hidden SU(2) gauge symmetry in our 
model. In the gapped B phase, the ground state has a non-trivial topological number, the spectral 
first Chern number or the chiral central charge, which reflects the chiral anomaly of the edge state. 
We proved that the topological number is identified to the reduced eta-invariant and this anomaly 
may be cancelled by a bulk Wess-Zumino term of SO(3) group through an index theorem in 2+1 
dimensions. 

PACS numbers: 75.10. Jm,03.67.Pp,71.10.Pm 



I. INTRODUCTION 

The concept of the topological order recently is widely 
interesting the condensed matter physicists because it 
may describe the different 'phases' without breaking any 
global continuous symmetry of the system However, 
unlike the conventional order related to the symmetry 
of the system in Landau's phase transition theory, the 
topological order of quantum states is not well defined 
yet. For example, in the quantum Hall effects, the topo- 
logical property of the quantum stats may be reflected 
by the filling factor of the Landau level which may be 
thought as a topological index, the first Chern number 
in magnetic Brilliouin zone @, H|. Nevertheless, only the 
first Chern number can not fully score the topological or- 
der of the quantum Hall states. In a given filling factor, 
the quasiparticlcs may obey cither abclian or non-abclian 
statistics. On the other hand, the edge state may par- 
tially image the topological properties of the bulk state 
In quantum Hall system, it was seen that the edge 
state may be described by a conformal field theory [5|. 
Thus, according to the bulk-edge correspondence due to 
the gauge invariance, it shows that the bulk state is de- 
termined by a Chern-Simons topological field theory Q. 
However, the bridge between the microscopic theory of 
the two-dimensional electron gas and the Chern-Simons 
theory was not spanned. 

Kitaev recently constructed an exactly soluble spin 
model in a honeycomb lattice Q- Using a Majorana 
fermion representation, he found the quantum state space 
is characterized by two different topological phases even 
there is not any global symmetry breaking. The A 
phase is a gapped phase which has a zero spectral Chern 
number and the vortex excitations obey abelian anyonic 



statistics. The B phase is gapless at special points of 
Brilliouin zone. When the B phase is gapped by a per- 
turbation, it is topologically non-trivial and has an odd- 
integer spectral Chern number. (We call the gapless B 
phase the Bl phase and gapped one the B2 phase.) Ki- 
taev showed that if the spectral Chern number is odd, 
there must be unpaired Majorana fermions and then the 
vortex excitations obey non- abelian statistics. Consis- 
tent with the non-abelian statistics, the fusion rules of 
the superselection sectors of Kitaev model are the same 
as those of the Ising model. However, the source of the 
non-abelian physics has not been clearly revealed yet. On 
the other hand, the first Chern number can only relate 
to an abelian group and therefore, an odd spectral Chern 
number leads to a non-abelian physics but an even one 
did not is topologically hard to be understood. 

Although Kitaev model has a very special spin cou- 
pling, its very attractive properties caused a bunch of re- 
cent studies % 1 EE EL Gl EI EI EE M EI El, El M ■ 

It is convenient to understand Kitaev model if one can 
map this model to a familiar model. In fact, Kitaev 
model may be mapped into a special p-wave paired BCS 
state if only the vortex-free sector of the model is con- 
sidered [ll| . We recently generalized Kitaev model to an 
exactly soluble model whose vortex-free part is equiva- 
lent to Ai x p x + AiyPy + i(A 2x Px + A 2y p y )-wave paired 
fermion states with tunable pairing order parameters A a b 
on a square lattice. [l2| . The phase diagram of our 
model has the same shape as that of Kitaev model, i.e, 
the boundary of the A-B phases are corresponding to the 
points p = (0, 0), (0, ±7r) and (±7r,0) in the first Bril- 
liouin zone. The A phase is gapped and may be iden- 
tified as the strong paring phase of the p-wave paired 
state [2l[. The B phase can be either gapped or gapless 
even if T-symmetry is broken. We find that gapless ex- 
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citations in the B phase, i.e., the Bl phase, is protected 
by a generalized inversion (G-inversion) symmetry under 

p x *-> -frf-Py and the emergence of a gapped B(B2) phase 
is thus tied to G-inversion symmetry breaking. For in- 
stance, the p x + ip y wave paired state is gapped while 
p y + ipy-w&ve paired state is gapless although they both 
break the T-symmetry. The critical states of the A-B 
phase transition remains gapless whether or not T- and 
G-inversion symmetries are broken, indicative of its topo- 
logical nature. Indeed, if all A a ;, are tuned to zero, the 
topological A-B phase transition is from a band insulator 
to a free Fermi gas. The Fermi surface shrinks to a point 
zero at criticality. 

In this paper, we further generalize the model pro- 
posed by the present author and Wang in ref. [l2| to 
a model whose square lattice mapping includes a next 
nearest hopping of the spinless fermions. In this case, 
the A phase is still a strong pairing phase as before. 
However, the B2 phase has more fruitful structure. The 
particle-hole symmetry is broken even if the chemical po- 
tential and the pairing parameters vanish. Near the long 
wave length limit (p* = (0,0) critical line), the effec- 
tive chemical potential has the different sign from that 
of the nearest neighbor hopping amplitude. Near other 
two critical lines (0, it) and (tt,0), when the next near- 
est neighbor amplitude is positive, the effective chemical 
potential is also positive. When the next nearest neigh- 
bor amplitude is negative, the effective chemical poten- 
tial is also negative. A positive chemical potential cor- 
responds to a closed Fermi surface of the particles and 
then a Pfaffian of the particles while a negative chemical 
potential to a closed Fermi surface of holes and then an 
anti-Pfaffian of the holes of the spinless fermion. There- 
fore, a Pfaffian/anti-Pfafhan phase transition happens in 
the B2 phase. This Pfaffian/anti-Pfafhan phase transi- 
tion has been seen in the context of the v = 5/2 fractional 
quantum Hall effect [22|, HH . The model we present here 
is exactly the same as a toy model on square lattice to 
study the Pfaffian and anti-Pfaffian physics [23]. The Bl 
and B2 phases when the next nearest neighbor hopping is 
absent are corresponding to the particle-hole symmetry 
is conserved or spontaneously broken. 

The another topic of this paper is trying to reveal the 
mathematical connotation behind the topological order. 
We emphasize that there is a hidden SU(2) gauge symme- 
try in this model if the model is represented by Majorana 
fermion operators. This non-abelian gauge symmetry is 
the source of the non-abelian physics of the model. The 
non-abelian degrees of freedom in the A phase are con- 
fined while in the B2 phase, the non-abelian degrees of 
freedom are deconfmcd. There is a Wess-Zumino(ZW) 
term for SU(2)/Z2 group whose lever k may character 
the confinement-deconhnement phases. A level k WZ 
term corresponds to a level k SU(2)/Z2 Chern-Simons 
topological field theory. It was known that k = 1 theory 
can only have abelian anyon while k = 2 theory includes 
non-abelian anyons [24] ]. A recently proved index theo- 
rem in 2+1 dimensions shows that the sum of this WZ 



term and a reduced eta- invariant fj is an integer [25[ . We 
show that difference between the WZ term and a part 
of the eta-invariant gives an ambiguity of the WZ term. 
Another part of this eta-invariant is identical to the chi- 
ral central charge, a half of the spectral Chern number. 
Thus, an odd Chern number corresponds to a 7ri ambi- 
guity while an even Chern number to a 27ri ambiguity. 
According to the bulk-edge correspondence, the former is 
consistent with k = 2 while the latter is consistent with 
k = 1. 

The rest of this work was organized as follows. In Sec. 
II, we recall Kitaev model and show the SU(2) gauge in- 
variancc. In Sec. Ill, we will describe the generalized 
model. In Sec. IV, we give the phase diagram of the 
system. In Sec. V, we consider the continuous limit of 
our model and show that the low energy effective theory 
is the Majorana fermions coupled to a SO (3) gauge field 
in a pure gauge. In Sec. VI, we apply the index theorem 
on odd manifold to our model. In Sec. VII, we present 
a understanding to the edge state from the index theo- 
rem point of view. The section VIII is our conclusions. 
We arrange three appendices. Appendix A is to address 
the mathematic expression of the index theorem on odd 
manifold because most of physicists are not familiar with 
it. In Appendix B, we give an introduction to the rep- 
resentation to the spin- 1/2 in the conventional fermion 
and Majorana fermion. And in Appendix C, for com- 
pleteness, we recall the vortex excitations in our model 
although it was studied in our previous work [T3 ] . 

II. KITAEV MODEL 

We first recall some basic results of Kitaev model, 
which is a spin system on a honeycomb lattice 0. The 
Hamiltonian is given by 

H ki = -j x e ofoj-jy E E 

x— links y — links 2— links 

where a a are Pauli matrices and 'x-,y-,z- links' are three 
different links starting from a site in even sublattice 0. 
This model is exactly solvable if one uses a Majorana 
fermion representation for spin. Kitaev has shown that 
his Hamiltonian has a Z2 gauge symmetry acting by a 
group element, e.g., for (123456) being a typical plaque 

W P = alalolololol 

with [Hki, Wp) = 0. In fact, we can show that this model 
has an SU(2) gauge symmetry in the Majorana fermion 
representation. Let & x ,y,z and c be four kinds of Majorana 
fermions with & 2 = & 2 , = bl = c 2 = 1 and define 

One observes SU(2) gauge invariant operators 

* a = ^Tr[x t x(a a ) T ] = l -(b a c - \e ahc b h b c ) (2) 
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cd 



with respect to the local gauge transformation X 
U cc ' x c ' d and then ( X ] ) cd -> ix ] ) cc ' {U' 1 ) ' 4 for 17 €SU(2) 
[26[. It is easy to check that a a /2 may serve as spin- 
1/2 operators. Replacing a a by er a , Kitaev model has 
a hidden SU(2) gauge symmetry which is trivial in the 
spin operator representation. The constraint D = 1 is 
also gauge invariant because D = b x b y b,c = —ia x a y a z . 
Under this constraint, a a takes the form ib a c after using 
b x b y b z c = 1. The SU(2) symmetry of er° = ib a c can be 
directly checked 



(-/i(p)M + A(p)<7* with A(-p) = -A(p), h = 
2 J + 4 J cos \p x cos and / 2 = 4 J cos \p x sin 

(/ = /i + if 2)- Assume ?/> (p) to be the solutions 
of Schrodingcr equation ifo(p 5 A)-0 ± (p) = ±E(p)ijj ± (p) 
with -B(p) = \/l/l 2 + |A| 2 . After normalization, we have 
L ■ aip = ±ip with 



V3JE(p) 



(-/a(p),-/i(p),A(p)). 



(7) 



where 



i6 a c = «6' a c', 



ai a 2 /?i ~/?2 

-CK2 ai — /?2 —A 

-/3i /3 2 «i a 2 

/3 2 /3i ~a 2 otx 




with 



a 2 + 0\ + /3| = 1 and b x b y b z c 



Kb'ybW 



(3) 



(4) 



1. 



Using Jordan- Wigncr transformation, a variety of Kitaev 
model on a brick-wall lattice has been exactly solved @ 
and a real space ground state wave function is explicitly 
shown [ll|. This variety should correspond to another 
gauge fixed theory. 

After some algebras, Kitaev transferred the Hamilto- 
nian to a free Majorana fermion one 0] 



H 



E 

p;fi,u— b,w 



H(p)auC- 



P,/X u p,I/5 



(5) 



Explicitly, near p = p* = -fpi + |p 2 mod(pi, p 2 ) 
with pi and p 2 the dual vectors of ni and n 2 , it is 
L = ^(6p^5p^^0) = (8p y ,Sq x ,A) with E(p) = 

A 2 /3J 2 . Near p = -p*, it is L = 
(6py,—6p x ,—A). According to Kitaev, one can define a 
spectrum Chern number by using the vector field L. We 
will be back to this issue when studying the index theo- 
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where H(—p) = —H(p) and c q /J are the Fourier com- 
ponents of a Majorana fermion operators and /Lt = b or 
w refers to the even or odd position in a z-link 0. The 
ground state is vortex-free and the corresponding Hamil- 
tonian Ho(p) is given by 



FIG. 1: Phase diag 
This is a (l,l,l)-cross section in all positive region. 



^o(p) 



i/(p) 
-if*(p) 



(6) 



with /(p) = 2( J x e lp ni + J y e l P' n2 + J z ). Here we still fol- 
low Kitaev and choose the basis of the translation group 
ni i2 = (±1/2, V3/2). In the next section, we will see 
that deforming the angle between the x-link and y-link 
to a rectangle will be much convenient. The eigenen- 
ergy may be obtained by diagonalizing the Hamiltonian, 
which is £b(p) = ±|/(p)|- The phase diagram of the 
model has been figured out in Fig. 1. Kitaev calls 
the gapped phase as A-phase and the gapless phase B 
phase. The A phase is topologically trivial and gapped. 
It is the strong-coupling limit of SU(2) like the antifcrro- 
magnetic Heisenberg model j26| and can be explained as 
the strong paring phase in the p wave sence [21[. After 
perturbed by an external magentic field, the B phase is 
gapped and has a non-zero spectral Chern number and 
then is topologically non-trivial citeki. Without lose of 
generality, wc consider J x = J y — J z = J. The effective 
Hamiltonian is then given by Ho(p, A) = (— f2{p))o' x + 



III. GENERALIZED EXACTLY SOLUBLE 
MODEL 




FIG. 2: (Color online): The honeycomb lattices and links. 



We now consider the Hamiltonian which is generaliza- 
tion of the Kitaev model in honeyconmb lattice to the 
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following one 

H = -J x £ a?a*-J y £ ofojf-J, E 

y— links 



this spin model and then the Hamiltonian reads 



x— links 
Kx ^2 <T b°b+e ]t <T b+e 3t +e lt 



Ex y z 

w 

Ez x y 
a b Pb+e^b+ez+ey 

b 

K V E/ (T w cr w+ey Cr w+ey+e z 



link* 



2^ <T b <T b+e z cr b+e z +e x a b+e z +e x +e z 



A\ ™ Z X X z 

y 2-^ ah a b+e z a b+e z +e y cr b+e z +ey+e z 
b 

- BbY^VbVb+eSb+ez+e^l+^+ev-ey 

b 

w / j <J w°-w+e :c (J w+e :c -ey (J w+e :c -e y -e z 
W 

- B w J2 <-e v °b°b+e z ^b+e z -e* 

b 

~ Bb E Vb-ey-e^b-eyVbVb+e^b+e^+e^b+e; 



(8) 



+e x +e z 



where '«/ and '6' labels the white and black sites in lat- 
tice and e x ,e y ,e z are the positive unit vectors, which 
are defined as, e.g., ei2 = e z ,e23 = e^egi = e y (See 
Fig. [2]). J x , y ,z, K x,y, ^x.y and Bb w are real parameters. 
This is a generalization of Kitaev model with the three- 
spin, four-spin and six-spin terms. It is easy to check this 
generalized Hamiltonian still has a Z 2 gauge symmetry 
acting by a group element, e.g., Wp = o\o\o\o\o\o% 
with [H 7 Wp] =0. In fact, one can add more Z 2 gauge 
invariant multi-spin terms, e.g., 

z V V V T. 

Z y V V Z V 

cr 9 cr 10 cr l cr 2 f7 3 cr 4 J 

z y v v v z 
cr 9 (T^cr^(T|(7|(7 16 , 

and so on. The site indices are shown in Fig. [5] For our 
purpose, however, we restrict on ([8]). 



FIG. 3: The effective square lattice. 
We now use the Majorana fermion representation for 



H 



— i E] Ja u ii C i C j ~ ^ El Kb,b+e z c b°b+e z +e 

a a— links b 

* / ] ^lo+e,- e x ,w — e z — e x c iu+e x — e x c io — e z — e x 
~ * E ^b.b+2e,+e 3 , C b C fe+2e 



; E A * 



2e z +e x 



-2e z — e x CwCw — 2e 



+ j/— partners 

+ » E /^6,6+e,+e.,-e 1 ,Ci,Ci, + e,+e a ,-e y 
6 

F ' &w,b+ey+e x c w c b+e y +e:c 

w 

+ * E Pt>,h+e,+ey+e m +e z CbCb+e z +e m +e v +e a (9) 
6 



Where ^ fe + ez = ^xW^+e^fe+e.+e^b+e^ A M+2e,+e x = 
^x U b,b+e z U b+e z ,b+e z +e w U b+e z +e n: ,b+e z +e w +e z 



and 



6+e J .' U '6+e 21 6+e 2 +e a , "fc+e = +e x: b+e 2 +e x -e y 



on a-links. 

y 



Pb, 



b+e z +e x —e v 



and 



t ~^w^w+e a; —ey—e z — 
M iD+e x -e H ,iu+e x -e y -e. °t C - 



B w U w : w+e x U w+e x ,w+e x — e y 

It can be shown that 
the Hamiltonian commutes with u% and thus the 
eigenvalues of Uij = ±1. Since the four spin and six-spin 
terms we introduced are related to the hopping between 
the 'b' and V sites, Lieb's theorem [27| is still applied. 
Following Kitaev, we take Ubw = ~u w f, = 1 and the 
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vortcx-frcc Hamiltonian is given by 

Hq iJ-x ^ ^ {CsbC-s — e x ,w Cs 7 wCs— e x ,b) 

s 



+ y partners + i J 2 c sb 

S 

s s 

5 S 

— iJx ^ ^ {CsbC-s — e x ,w Cs,wCs—e x ,b) 
s 

~t~ ^ ^ {Cs,bCs — e x ,w H~ Cs,wC-s—e x ,b) 

s 

~t~ ^~0~ ^ ^j ^s,b^s-\-e x ,b H~ Cs.-iD Cs— ,tt; ) 
s 

+ y partners + i J z ^ c 5& c sw 

s 

^ ^ (^s.&^s+e^ — ,w Cs,wCs-\-e x —e y ,b) (^^) 

2.6 ^ ^ (^s^Cs-fe^ — e y ,iu "I - ^s.^Cs-j-e^ — e y ,b) 
s 

~l~ ^^( < Cs,bCs J re x -\-e y ,w Cs.ttjCs+e^ +e y ,&) 

3 



where s represents the position of a z-link, A a 
and J a = 



J Q +A 



J a -X a g± _ B b ±B u 



2 " 2 

To simplify the pairing, one takes £? + — and denotes 
B = B~ . Defining a fermion on z-links by [ill, 

d s = (c s ,b + ic s ^ w )/2, d\ = (c Sjb - ic s ^ w )/2, (11) 

the vortex-free Hamiltonian becomes an effective model 
of spinless fermions on a square lattice (Fig{3]) 

Ho = J.^(dK-l/2) 

S 
S 

+ J x (d\d s+ex - d s d\ +e J 

+ \ x J2( d l+ejl - d s+ejs) 

S 

+ iK x ^2{d s d s+ex +S s d\ +e J + y partncrs.(12) 



Or it is 

H = £(-*dJrfj + + h.c.) - ^(4% - 1/2) 

(ij) 1 

-t (4di+^c)+s52(4d i ± em -4d i ± ey ), (is) 

({ij}) i,± 



where t = — J "t v , t' = —B, /j, = —J z and S 
The paring parameters are defined by Ai t i± en . = \ x ^ y + 
in x ^ y . The last equality in eq. (fT3]) is the toy model Hamil- 
tonian describing Pfaffian/anti-Pfaffian states [22]. Note 
that the pairing free Hamiltonian is particle-hole sym- 
metry if Aij ~ /i = t' = 0. The t'-tcrm breaks the 
particle- hole symmetry. The <5-term breaks the it/2 ro- 
tational symmetry. After Fourier transformation d s = 
T^I-Ep e lp s rf P , we have 



Ho = J2^ d l d P + ^( d l d - P + d P d -p) 
p 



i^(dldl p -dpd_ p ) 



(14) 



2 v P - P 
where the dispersion relation is 

£,p = Jz — Jx cosp x — J y cospy + 2B cos p x cospy (15) 
and the pairing functions arc 

Ap = Ai a sinpz + Ai v sinp y , 



A p = A 2x sinpz + A 2y sinp y 



(16) 



with A 1)X(j/ ) = K x ( y ) and A 2 . x ( y ) = K( y )- After Bogoli- 
ubov transformation, 



a p = Updp - v p d_ p , 
«p = u * P d l - v P d -P 



the Hamiltonian can be diagonalizcd 

Ho = i? p apa p + const. 



(17) 



(18) 



and the Bogoliubov quasiparticlcs have the dispersion 

E p = v /^ + (Ap)2 + (A p )2. (19) 

The Bogoliubov-de Gennes equations are given by 

EpU p = ( p u p - A p v p , EpVp = -£p?> p - Apitp (20) 



with 



Vp/up = -(E p -£ P )/Ap, 
I"pI 2 -^(1 + ^), 

^ 2 = ^-wJ- 



(21) 
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IV. PHASE DIAGRAM 

We now study the phase diagram in parameter space. 
The phase diagram when t' = has been discussed 
in our previous work [l2| . which has the same shape 
as that in original Kitaev model (with (J Xl J y , J z ) in 
Fig.l substituted by {J x ,J y ,J z j) but the structures 
of the B phase are more fruitful. After including 
the i'-term, the phase boundary is still in (p Xl p y ) = 
(0, 0), (0, ±7r), (±tt, 0), (±7r, ±7r) as we know before. For 
the present model, it is J z ±J x ±J y = with J z = J Z +2B 
for (0,0),(±7r,±7r) and J z - 2B for (0, ±tt), (±tt, 0). 
In {J x ,J y ,J z ) space, the phase diagram are of the 
same shape as that of original Kitaev model (Fig. [TJ 
(Jx, Jy, Jz) is replaced by (J x , J y , J z ) ). The A phase 
is a strong pairing phase. The nature of the B phase 
is much more intriguing. Inside the B phase, £ p , Ai. p 
and A 2jP can be zero individually. The gapless condi- 
tion (E p — 0) requires all three to be zero at a common 
p*. This can only be achieved if (i) one of the A ap = 
or (ii) Ai p oc A 2 . p . If either (i) or (ii) is true, £ p and 
A p can vanish simultaneously, i.e. E p = at p*, and the 
paired state is gapless. Otherwise, the B phase is gapped. 
Note that contrary to conventional wisdom, T-symmetry 
breaking alone does not guarantee a gap opening in the 
B phase. The symmetry reason behind the gapless con- 
dition of the B phase becomes clear in the continuum 
limit where E p = implies that the vortex-free Hamil- 
tonian must be invariant, up to a constant, under the 
transformation p x «-> rjp y and J x «-> rf~ 2 J y with r/ = ^ a,y 
with a = 1 or 2 and for nonzero A. We refer to this 
as a generalized inversion (G-inversion) symmetry since 
it reduces to the usual mirror reflection when r\ = 1. 
This (projective) symmetry protects the gapless nature 
of fermionic excitations and may be associated with the 
underlying quantum order [28j . Kitaev's original model 
has Ai.i = 0, and is thus G-inversion invariant and gap- 
less. The magnetic field perturbation 0] breaks this G- 
inversion symmetry and the fermionic excitation becomes 
gapped. 

The continuum limit takes place near the critical line 
(0,0). In this case, 

J z -J x -J y - 2t' = -p eff - 2t' = 0. (22) 

with Me// = Jz — Jx — Jy Slight inside of the B phase, 

Me// > -2/'. (23) 

If t' < 0, Me/ / is positive and the system is in the Pfaffian 
phase. If t' > 0, 

Me// < 2\t'\ (24) 

may be either positive or negative. Me// > means that 
the electron Fermi surface is closed (Ay = 0) while the 
hole Fermi surface is opened. This is the d-particlc paired 
phase. On the other hand, /i e // < means that the elec- 
tron Fermi surface is opened while the hole Fermi surface 



is closed. This is the d-hole paired phase. Thus, there is 
a Pfaffian/anti-Pfaffian transition as Me// is across zero. 

There are other two critical lines (0, ir) and (it, 0) near 
which there are also gapless excitations. The critical con- 
dition is given by 

J z T Jx ± Jy - 2B = T $ + 2t' 

= 2f'-Me// = 0. (25) 

Inside of the B phase but near the critical lines, 

Me// - 2f > 0. (26) 

To satisfy this condition, Me// has to be the same sign 
as that of t'. Therefore, if t' > 0, the system is in the 
Pfaffian state and the system is in anti-Pfaffian state if 
t' < 0. 

We see that when t' < 0, the sign of the effective chem- 
ical near the critical line (0, 0) has the opposite depen- 
dence on the sign of t' to the effective chemical near other 
two critical lines. For t' > 0, the sign of the effective 
chemical near the critical line (0, 0) may change from the 
opposite to the same as that near other two critical lines. 
Hence, if t' ^ 0, there must be a Pfaffian/anti-Pfaffian 
transition inside of the gapped B phase. 

At t' = 0, the Pfaffian and anti-Pfaffian states are 
degenerate. As we have seen [l2[, in the gapless B 
phase, there are two gapless Majorana excitations at 
nodal points while in the gapped B phase, the particle- 
hole symmetry is spontaneously broken, i.e., the ground 
state is either Pfaffian or anti-Pfaffian. All above discus- 
sions are consistent with those in ref. 12211. 



V. CONTINUOUS LIMIT, DIRAC EQUATIONS 
AND SO (3) GAUGE THEORY 

In fact, the ground state wave function for a general p- 
wave paired state can also be calculated in the continuous 
limit. The BCS wave function is given by 

in) = ni u pl 1/2cx p(^E-9p rf P d - P )l°)' (27) 

p p 

where g p = v p /u p . For even fcrmion number N, the 
Pfaffian ground state wave function reads 

N/2 

*(ri,...,n)«5^BgnPjJff(r ftt _ 1 -r flM ) (28) 

P i=l 

with g p is the Fourier transform of g(r). For the A phase, 
g p oc A p and the analyticity of <7 p leads to g(r) oc e~ Mr 
as the same calculation in a pure p x + ip y strong pairing 
state. In the B phase, if the G-chiral symmetry is broken, 
defining p' a — A a bPb with a = 1,2 and b = x,y, g p oc 
jr^r and then g(r) = ^-^r with x' a = A^x b . This is 

a Pfaffian state with z' = x\ + ix' 2 and is corresponding 
to a weak paired gapped fermion state [2l| . 
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In the long wave length limit (small p limit), we ap- 
proximate £ p f» — m = J z — J x — J y and define ip{t, s) = 
(u(t,s),v(t,s)) with v) the Fourier transformation of 
(itp, Up). The BdG equations in the gapped B phase reads 

e^lad^ + irmp = (29) 

with the vielbein e M * = S^t and = Ay. 70 = a 2 
and 71 = a z a l . Furthermore, it may be generalized to 
a curve space with spin connection term added [2ll ]. In 
general, the equation of w is not compatible to that of 
v and the fermions are Dirac fermions. However, it is 
easy to show that for this p-wave paired state, u and 
v* obey the same equation, i.e., the anti-particle of the 
quasiparticlc is itself. Then the fermions are Majorana 
ones. 

We now turn to local gauge transformation. As we 
have pointed out, taking = 1 means the whole SU(2) 
gauge symmetry is fixed. If making an SU(2) gauge 
transformation, we will run out the vortex-free state. 
Keeping in the vortex-free state, the gauge transforma- 
tion is confined in an SU(2)/Z2 ~ SO(3) one. Therefore, 
making an SO (3) transformation, 

tf(T,t) = U(T,t)1>(T,t), (30) 

Dirac equations read 

e^ a D^' + imi* 1 = 0, (31) 

where the covariant derivative Dp = dp + Ap with Ap = 
U~ 1 dpU with U GSO(3). This gauge potential is a pure 
gauge with a vanishing strength F = dA + A A A = 0. 

We now discuss the nature of the gapless B phase 
in the general model with G-inversion symmetry. In 
this case, E p = at p = ±p* which are the solu- 
tions of £ p = and, say, A p = Ai p = 0. At p*, 
the fermion dispersions are generally given by 2D Dirac 
cones. However, by a continuous variation of the pa- 
rameters, one can realize a dimensional reduction near 
the phase boundary where the effective theory is in fact 
a 1+1 dimensional conformal field theory in the long 
wave length limit. Let us consider parameters that are 
close to the critical line with |sinp*| <C |cosp*| where 
,9 q = sgn[g x Ai x cosp* + q y A ly cosp*] = sgn(g£) with 
q = p — p*. Doing the Fourier transform, we find 

g{v) = J d q ' x d^ x ' + <y^ x ) 

= %') [dq' x -^ S mq> x x>~ 6 -^. (32) 

The <5(y)-function indicates that the pairing in the gap- 
less B phase has a one-dimensional character and the 
ground state is a one-dimensional Moore-Read Pfaffian. 
The BdG equations reduces to 

id t u = —iAi x (l + irj)d x iv, id t v = iAi x (l - ir))d x iu{33) 

with rj = ■^ JL . Thus, the gapless Bogoliubov quasipar- 
ticles are one-dimensional Majorana fermions. The long 



wave length effective theory for the gapless B phase near 
the phase boundary is therefore the massless Majorana 
fermion theory in 1+1-dimcnsional space-time, i.e. a 
c = 1/2 conformal field theory or equivalently a two- 
dimensional Ising model. 

VI. TOPOLOGICAL INVARIANTS AND INDEX 
THEOREM 

A. Spectral Chern Number and r/-invariant 

We note that there is no spontaneous breaking of a 
continuous symmetry in the phase transition from A to 
B phases. Kitaev has shown that the A phase in his 
model is topological trivial and has zero spectral Chern 
number while the gapped B phase has this Chern number 
±1 0] . This fact was also already realized by Read and 
Green in the p x + ip y paired state. Here we follow Read 
and Green [2l| to study this topological invariant for a 
general p-wave state. 

In continuous limit, p = (p x ,Py) is in an Euclidean 
space R 2 . However, there is a constraint |up| 2 + |w p | 2 = 1, 
which parameterizes a sphere S* 2 . As |p| — -> 00, £ p — > E p . 
Then, v p — > as |p| — > 00. Therefore, we can compact 
R 2 as an S 2 by adding 00 in which v p — * to R 2 . The 
sphere |u p | 2 + \v p \ 2 = 1 can also be parameterized by 

n p = (A { p \-A p 2 \£ p )/E p because |n p | = 1. (u p ,v p ) 
describes a mapping from S* 2 (p € R 2 ) to S 2 (spinor 
I rip I = 1). The winding number of the mapping is a 
topological invariant. The north pole is Up = 1,V P = at 
|p| = 00 and south pole u p = 0, v p = 1 at p = 0. For n p 

parametrization, n = (0, 0, ^) at |p| = 00 and (0, 0, ^jr) 
at p = 0, either the north pole or south pole. 

For strong pairing phase, we know that u p — * 1 and 
v p — > as p — > (or equivalently, £ p > 0). This 
means that for arbitrary p, (u p , v p ) maps p-sphere to the 
up-hemishperc and then winding number is zero. That 
is, the topological number is zero in the strong pairing 
phase. 

For the weak paring phase, u p — > and v p — * 1 as 
p — > 0. This means that the winding number is not zero 
(at least wrapping once). For our case, it may be directly 
calculated and the winding number is given by 

v = J~ J dp x dp y d Px n p x d Py n p • n p = 1 (34) 

Defining -P(p) = ^(1 + n p • a), which is the Fourier com- 
ponent of the project operator to the negative spectral 
space of the Hamiltonian, this winding number is iden- 
tified as the spectral Chern number defined by Kitaev 


v =^~j Tr(P-(d Px P-d Py P-d Py P_d Px P_)]dp x dp y (35) 

There is a topolo gica l invariant called Atiyah-Padoti- 
Singer eta-invraiant [33j which reflects the asymmetry of 
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the spectrum of the Dirac operator 
1 



fj(S 2 



lim 

2 s^o 



dXp(X)sgn(X)\X\- 



(36) 



where p(A) is the spectral density. Transforming the vari- 
able from A to p, the measure of the integration from A p 
to q includes a Jacobian determinant [29[ , 



4tt 



= i-/dV(A p ,p)i-. 



U [ d 2 pTr[J(n p -a,q)...} 
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Note that — ^n p • a is the signature matrix sgn(i/(p))of 
the Hamiltonian, one has 



fj(S 2 ) = ^Vimj dAp(A)sgn(A)|A|- 



dpTr[isgn(i7(p))p(isgn(i/(p)))],(37) 



where 



P(^(q))) - ^fM^ftpW 



9sgn(JJ(q)) 9sgn(#(q)) 



dq y 



dq x 



(38) 



Defining the spectral projector I — P = P-(q) = |(1 — 
sgniJ(q)), the eta-invariant is exactly equal to one half 
of the spectral Chcrn number defined by Kitaev Q 

V(S 2 ) = ^-J Tr[P_(q)dP_(q) AdP_(q)]. (39) 

Kitaev identifies one half of i/ as the chiral central charge 
c_. Our result shows that this chiral central charge is 
just the eta- invariant. Physically, it is easy to be under- 
stood because both c_ and fj reflect the anomaly of the 
spectrum of the system. 



B. Index Theorem 

In the continuous limit, if the space is compacted as 
S 2 , the 2+1 space-time is a ball X with a boundary B = 
S\{t = 1) U Y U S 2 {t = 0) where 5| are the top and 
bottom halves of a sphere and Y is a cylinder. Now, 
we can apply the index theorem (|A2ll to this spectral 
problem of the Dirac operator in X with boundary B 
|25| |. The general form of thye index theorem in an odd- 
dimensional manifold is briefly reviewed in Appendix A. 
In 2+1-dimensions, the index theorem for the Toeplitz 
operator reads 

IndT ff = -L f T r [(gdg- 1 ) 3 ]-f 1 (B,g) + T,(B,P,c^Q) 



where the first term is equal to T/2ni with T the WZ 
term. The Maslov triple index r fl (B,P,g) is an integer 
[39l ] . We do not have a physical explanation of r M (B, P, g) 
yet and it possibly relates to the central charge of the 
theory [40| . fj is the reduced eta-invariant. The first term 
in ([40|) determines the bulk state topological properties 
and the latter two terms reflect the boundary topological 
properties. That the index Ind T g — is an integer 
determines the bulk-boundary correspondence. 

The WZ term T is defined for the fundamental rep- 
resentation of SU(2) but g is restricted to a subgroup 
SU(2)/Z2 ~ SO (3). The reduced eta- invariant is given 
by 

f}(B,g) = St=i) +v(S-,9t=o) 



[fj(Sl,g T=1 )-fj(S 2 + ,9 T = )] 
+ [v(S+,9t=q) +fj(S 2 _,g T=0 )] 
Afj(S 2 + )+fj(S 2 ,g T=0 ), 



(41) 



because ij(Y) = for Y may contract to a cycle S 1 . 
In general, fj(S±,gi,o) =|[dim(ker£)(5|,5i i o)) + 
77(£)(5±,g , i i o))]- Because of a non-zero gap, 
dim(ker£>(5|)) = 0. Therefore, fj(S 2 ,g) = 
±Ar)(Sl) + \-q(D{S 2 ,g T=0 )). The discrete eigen- 
states of the Dirac operator do not contribute to fj 
because there is no asymmetry of the spectrum for these 
states. fj(S 2 ,g T= o) = fj{S 2 ) is just the eta-invariant 
calculated in the previous subsection. 
Now, the index theorem reads 

Ind T g = r/2ni - v/2 - Afj(S 2 + ) + t^B, P, g). (42) 

The integrity of Ind T g — t^(B, P,g) requires 

AT/2ni ee T/2Tri - Afj(Sl) = v/2 mod(Z). (43) 

Dai and Zhang have thought Afj(S+) as an intrinsic form 
of the WZ term [25| and then Ar is in fact an ambiguity 
of the WZ term . It is seen that, if v is even, one requires 
Ar = 27rjx integer. If v is odd , AY j2-Ki is required 
to be a half integer. In our model, it is known that for 
an SO(3) group, Ar = 7rixinteger, which is consistent 
with v = 1. In general, the index theorem (|40|) gives a 
constraint to the WZ term. An odd v requires an even 
level k WZ term in the effective action and the minimal 
one is k = 2. This is consistent with the non-abclian 
anyonic statistics of the vortex excitations. For an even 
v, the minimal value of k is one, it is consistent with the 
abelian anyonic statistics. 



VII. EDGE EXCITATIONS AND BULK-EDGE 
CORRESPONDENCE 

In the previous discussion, the two-dimensional space 
is taken to be a tours without boundary. If we consider 
a two-dimensional space with edge instead of the torus, 
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Kitaev has shown that the gapless chiral edge excita- 
tions coexist with a non-zero spectral Chern number 
This is a general result if the bulk states are gapped. 
In fact, it was generally known that the eta-invariant of 
the Dirac operator can be related to the ground state 
fermion charge 0, [53, OHl- By using the continuous 
equation dtp + V • j = 0, the eat-invariant is related to 
the edge current integrated along the one-dimensional 
edge S, i.e.,[I| 

r, = 2 [ ds j , (44) 
Js 

where the factor '2' is different from T' in (104) of Kitaev 
in Q because of \ factor in ([5]). Thus, a non-zero eta- 
invariant corresponds to a non-vanishing net edge cur- 
rent and then the gapless chiral edge excitations. Hence, 
the index theorem already explained the bulk-edge cor- 
respondence. 



VIII. CONCLUSIONS 

We studied a generalized Kitaev model whose vortex- 
free sector can be mapped to a p-wave paired state with 
the next nearest neighbor hopping. The phase diagram 
is figured out. The property of the gapped B phase was 
very interesting for a Pfafuan/anti-Pfaffian phase tran- 
sition was found in this phase. According to the gauge 
invariance of the spin- 1/2 theory in the fermion repre- 
sentation, we found the low-lying effective theory of the 
model is described by Majorana fermion coupled to a 
gauge field. The existence of this non-dynamic gauge 
field enabled us to understand the mathematic connota- 
tion behind these topological orders. The edge conformal 
anomaly can be cancelled by the bulk WZ term in terms 
of the recently proved index theorem on odd manifold. 
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APPENDIX A: INDEX THEOREMS: 
MATHEMATICAL PREPARATION 

Index theorems relate the analytical index of a differ- 
ential operator to the topological index of a vector bundle 



that the operator acts on. Atiyah and Singer (AS) proved 
a general form of the index theorem for even dimensional 
compact manifold [H)]. It said that the analytical index 
of an elliptical differential operator on the vector bundle 
is a global topological invariant which can be expressed 
by the integral of local topological characters on the back- 
ground manifold. A generalization of the theorem to even 
manifolds with boundary, so-called Atiyah-Patodi-Singcr 
(APS) theorem [33|, e.g., for a Dirac operator D on a 
manifold M with boundary dM , reads 

IncLD = / A(M) + -(h[dM\ + rj[dM]) + u(dM),(Al) 

where A is the Hirzebruch ^4-class of M, h is the di- 
mensions of the zero modes of the boundary Dirac op- 
erator Dqm and r\ is the APS eta-invariant defined by 
7](dAI) = lim s ^ Ea^o s g n MI^I _s where A is the eigen- 
value of Dqm ■ lo{8M) is a Chcrn-Simons term caused by 
the non-product boundary metric and physically giving 
anomaly in quantum field theory [34| . There are many 
applications of the AS and APS theorem in physics, e.g, 
see Refs.p9l I30L l3l|. Recently, index theorems were ap- 
plifid to chiral p-wavc superconductors [35| and graphene 

There are partners of AS and APS theorems on odd- 
dimensional manifolds X (tf=odd)[25|,[37],[38j|. Again, we 
consider Dirac operator D. C = J2\ is the spectrum 
space of D where C\ is the subspace with eigenvalue A. 
P is a project operator defined by PL = C+ = X)a>o 
Let C trivially take its value on C , i.e., a state if,' E C 
is extended to an TV vector which transforms under a 
group GL(N,C). For Toeplitz operator T g = PgP on 
odd-dimensional manifolds X with boundary B in which 
the group element g is not an identity, an index theorem 
is given by [25[ 

Ind T g = - 1 f A(R x )Tr[cM-R C )]Mg) 
(2m) 2 Jx 

- f][B,g]+T IJl (B,P,g), (A2) 

where g £GL(N, C) (or a subgroup like SU(A r )); R x and 
R c are the curvatures of the background manifold X and 
£; fj [B,g] is a reduced eta-invariant and r^, which is an 
integer, is the Maslov triple index (39[. We do not have 
a physical explanation of r M yet and it possibly relates to 
the central charge of the theory [4(J . ch(g) is odd Chern 
character defined by 

ch (») = E (2^T^ Tl[{g ~ ld9)2n+1] - (A3) 

n=0 ^ '■ 

The first term in (|A2[) determines the bulk state topologi- 
cal properties and the latter two terms reflect the bound- 
ary topological properties. That the index Ind T g — r M 
is an integer determines the bulk-boundary correspon- 
dence. 
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APPENDIX B: SPIN AND MAJORANA 
FERMIONS 

In this appendix, we introduce the Majorana represen- 
tation of spin-1/2 operators. Consider the Pauli matrices 



1 

1 



,<7 a 



-i 

1 



, a 



1 

-1 



. (Bl) 



with a x a y — ia z and so on. The spin-1/2 matrices are 
one-half of the Pauli matrices: S = §. Casimir of SU(2) 
group requires S ■ S = S(S + 1) = 3/4 , i.e., a ■ a = 3. 
Using the conventional fermion operators cj and cj., the 
spin-1/2 operators can be expressed by 



That is, 

b x = cj +c T , W = i(c\ -c T ), 
b z =c\ +C|, c = i(c\ -c±). 

It is easy to check that 

b 2 a = l, c 2 = l 

bah = -b b b a ,cb a = -b a c 



(B6) 



Therefore, using the Majorana fcrmions, one can express 
the spin operators by 



a x = c\c d +cjc T 


i 

a x = ^{b x c-b y b z , 


a y = -i(c\c l - cjc T ) 


av= l -(b y c-b z b x ), 


= 4 C T - C W (B2) 


v z = -j{b z c - b x b y ). 


That is a a = 4<j%,c s >. Since {c,^} = l,c 2 = c^ 2 = 0, 
one has a x a y = icr 2 . However, it is easy to see that 


Correspondingly, the constraint reads 


(<T a ) 2 = n T - 2n- l n l 


T X = l ~{b X C+byb Z ) = 


. To insure the Casimir operator constraint, one requires 


T V = l -{byC+b z b X ) = 


n T + nj, = 1, ntnj, = 
These are equivalent to an SU(2) constraint 


T z = l -{b z c + b x b v ) = 



(B7) 



T* = 4cJ + Ci c T =0, 
T y = i(c[cj — c;C|) = 0, 
T z = n T + n L - 1 = 



(B3) 



with T x T y — iT z and so on. It is well-known taht 
{& a /2,T a /2} form anSO(4) ~ £17(2) x SU(2)/Z 2 Lie 
algebra. 

The fermion expression brings extra degrees of freedom 
and then there is an SU(2) gauge invariant of the spin 
operators. Defining 



(Xa/3) 



the spin operator may be rewritten as 



1 



* a = 2 Tr [x f x(^) T ] 



(B4) 



(B5) 



Making a gauge transformation Xap ~ * 9afX-y/3 an d 

then xlp XaW-rP for 9 eSU ( 2 ) and .9 f = .9~\ onc nas 
<r a is gauge invariant. 

Majorana fermions arc related to the conventional 
fermion through 

c T = -^(bx-iby), c L = hb z -ic) 



(B8) 

Compactly, it is D = b x b y b z c = 1. In an explicit gauge in- 
variant form, D = —i<j x o v <j z = 1. Using this constraint, 
the spin operators are simplified as 

a a = ib a c. (B9) 



APPENDIX C: ISOLATED VORTEX 
EXCITATIONS 

The isolated vortex excitations above the Pfaffian 
ground state has been discussed in [u|. For self- 
containing of the paper, we repeat the paragraph which 
discuss the vortex excitations. The Z 2 vortex excitation 
in the spin model which corresponds to setting Wp = — 1 
for a given plaquette. Although the Hamiltonian in the 
fermion representation is bilinear, it is difficult to obtain 
analytical solutions of the wave function with vortex ex- 
citations . Our strategy is to evaluate the energy of the 
Moore-Read Q trial wave function with two well sepa- 
rated half-vortices located at wi and w 2 shown in Fig. 4 
. In the gapped B phase with p x + ip y -w&ve pairing, 

*(zi, ...z N ;wi,w 2 ) oc Pf(g'(z i ,Zj;w 1 ,w 2 )), 

(zi - wi)(z 2 - w 2 ) + (wi «-» w 2 ) 



g'{zr,z 2 )wx,w 2 ) cx 



2 V " yn 1 2 
\=?Q> x + ib y ), c\ = -(b z +ic) 



z\ - z 2 



\w!,w 2 ) tx exp{i ^ 9\zx,Z2\w\,w 2 )d\ 1 d\ 2 } (Cl) 
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Performing a Fourier transformation, we have 

\ Wl ,w 2 ) ocexp{^ fl UK)4 +k 4_ k }, (C2) 



K,k 



where k = ki — k 2 and K = ki + k2 are the relative and 
the total momenta of the pairs and <7j.(K) is the Fourier 
transform of g'(rx,T 2 ). One can show that, in a system 
with linear dimension L, g' k (K = 0) ~ r(l/6 + i/8 — 
(1 + i)(wi + w 2 )/8L + W1W2/L 2 ). The Hamiltonian in 
the presence of the two vortices shown in Fig. 1 where 
the red z-links have Ub w = —1 and all others Ub w = 
1, may be written as H = Hq + SH. Here Hq is the 
vortex-free Hamiltonian and SH is the vortex part. The 
latter is expressed as a sum of (twice) the pairing and 
chemical potential terms in Eq. (6) over the red z-links 
extending in the ^-direction (the line with x — y) between 
the vortices. By a direct calculation, one can prove that 
SH has the following form 



(w 1 ,w 2 \SH\w 1 ,w 2 ) 
i( e »»i(Ps+Pe) 

— 



(C3) 



-f(Pi,^) = o, 



where f(p^,p'^) is an analytical function of p^+p't 
other hand, one can check that since [Hq,J2k kd'ki^- 



0)<%+k^k-k] = K ^ sector does not play a 

nontrivial role in calculating E v . The energy of such a 
vortex pair is given by 

K = (w 1 ,w 2 \H\w 1 ,w 2 ) = (w 1 ,w 2 \H \w l ,W2) 
= ^2 E k \u k 5g k \ 2 (wiw 2 \d k dl}w 1 W2) 



(C4) 



where <5c/k = g k (K = 0) — g^ and A(w±, w 2 ) is a positive 
constant. In general, (wiw 2 \dkd k \wi i w 2 ) = 1 — |g k (A" = 
0)| 2 is the quasiholc distribution when the two vortices 
are located at wi and w 2 . Thus, E v is indeed the en- 
ergy cost to excite the vortex pair. Since g k ~ 1/k and 
Ek ~ k for fc — * 0, there is no infrared divergence in E v . 
Therefore it costs a finite energy for such vortex pair ex- 
citations which is isolated to other states in the gapped 
B phase. Note that A(wi,w 2 ) — > as wi and w 2 — > 00, 
which corresponds to the case where the bulk of the sys- 
tem is sent back to the vortex-free ground state. 




FIG. 4: The vortex excitations. The grey solid circles denote 
Wp — 1 and the red solid circles denote vortices with Wp = 
-1. 

The hniteness of E v and vanishing of (SH) = im- 
ply that the vortex excitations are isolate either to the 
ground state or the other excitations. This is consistent 
with analysis of Read and Green to the vortex excitations 
in U(l) vortex excitation of the p-wave paired state [21 [. 
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